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The Sunyaev-Zeldovich (SZ) effect was previously measured in the Coma 
cluster by the Owens Valley Radio Observatory and Millimeter and IR Testa 
Grigia Observatory experiments and recently also with the Wilkinson Microwave 
Anisotropy Probe satellite. We assess the consistency of these results and their 
implications on the feasibility of high-frequency SZ work with ground-based tele- 
scopes. The unique data set from the combined measurements at six frequency 
bands is jointly analyzed, resulting in a best-fit value for the Thomson optical 
depth at the cluster center, r = (5.35 ± 0.67) x 10~ 3 . The combined X-ray and 
SZ determined properties of the gas are used to determine the Hubble constant. 
For isothermal gas with a (3 density profile we derive H = 84 ±26 km/ (s ■ Mpc); 
the (la) error includes only observational SZ and X-ray uncertainties. 

Subject headings: cosmology: cosmic microwave background - observations - 
galaxies: clusters: individual (A1656) 
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1. Introduction 

The Sunyaev-Zeldovich (SZ) effect (Sunyaev & Zel'dovich 1972) constitutes a unique 
and powerful cosmological tool (for reviews, see Rephaeli 1995a, Birkinshaw 1999, Carlstrom, 
Holder, & Reese 2002). Many tens of high quality images of the effect have already been 
obtained with interferometric arrays operating at low frequencies on the Rayleigh- Jeans side 
of the spectrum. Multifrequency SZ measurements were made of only very few clusters, so 
the potential power of spectral diagnostics has not yet been sufficiently exploited to reduce 
signal confusion and other errors. Systematic uncertainties are the main hindrance to the 
use of the effect as a precise cosmological probe. These can be most optimally reduced when 
high spatial resolution X-ray and spectral SZ measurements of a large sample of nearby 
clusters are available in the near future. The multifrequency capability of many upcoming 
ground-based bolometric array projects is crucially important for the separation of the SZ 
signal from the large confusing atmospheric signals. 

The SZ effect was measured in the Coma Cluster with the ground-based Owens Val- 
ley Radio Observatory (OVRO; Herbig, Lawrence, & Readhead 1995) and Millimeter and 
IR Testa Grigia Observatory (MITO; De Petris et al. 2002) telescopes. The Wilkinson Mi- 
crowave Anisotropy Probe ( WMAP) team has recently reported measurement of the effect in 
Coma at 61 and 94 GHz (Bennett et al. 2003b); this is the first measurement of the SZ effect 
from a satellite. Together, these measurements yield the first SZ spectrum with six spectral 
bands. The relatively wide spectral coverage obtained with ground- and space-based tele- 
scopes allows a much-needed gauge of the impact of atmospheric emission on high-frequency 
ground-based observations of the SZ effect. The quality of the spectral results allows a 
meaningful determination of the Hubble constant, H . Even though this important parame- 
ter is thought to have been determined quite precisely from the WMAP first-year sky survey 
(Bennett et al. 2003a), it is, of course, very much of interest to measure it independently by 
various methods. The determination of H from cluster SZ and X-ray measurements yields 
additional insight beyond what is offered by an all-sky cosmic microwave background (CMB) 
anisotropy survey: in principle, the cluster SZ & X-ray (SZ-X) method provides a test of 
both the isotropy of the cosmological expansion and its temporal character. In this Letter we 
briefly discuss the qualitative deductions that can be drawn from the general agreement be- 
tween the ground-based and WMAP results, and use the full database on Coma to determine 
H . 
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2. SZ Measurements of the Coma Cluster 

The rich nearby Coma Cluster has long been a prime target of SZ observations (e.g., Birkinshaw, 
Gull, & Northover 1981, Silverberg 1997); the effect was detected towards Coma with the 5.5 
m OVRO telescope, with a 7'. 3 beam, operating at 32 GHz. A long series of drift scans with 
the MITO 2.6 m telescope, which had a 16' beam, led to the detection of the effect at three 
spectral bands centered on 143, 214, and 272 GHz. Even though the WMAP telescope is not 
optimized for SZ observations, the effect was detected with the W- and V-band radiometers 
at 61 and 94 GHz, with beam sizes of 20' and 13', respectively (Bennett et al. 2003b). The 
predicted SZ size of Coma is more than ~ 30', so these measurements were not substantially 
affected by beam dilution. Results of the measurements are listed in Table 1; these reflect 
a slight revision of our previously reported SZ signals as result of an improved calibration 
of the the MITO measurements (Savini et al. 2003). The SZ spectrum of Coma is shown in 
figure 1. 

Note that the MITO error bars are about a factor 3 smaller than those of the WMAP 
data; this is mainly due to the difference in observing time. Given the low signal-to-noise 
ratios (S/Ns) of the two WMAP data, the main advantage of these satellite observations lies 
in the capability of relatively precise calibration through measurements of the modulation 
of the CMB dipole anisotropy. This is estimated to be ~ 0.5% by the WMAP team, as 
compared with ~ 10% for MITO and from ~ 3% to ~ 10% for OVRO. It is interesting 
that the relative calibration of MITO and WMAP are in agreement to within ~ 10%. This 
result provides further support for ground-based SZ observations, which are susceptible of 
substantial uncertainties because of atmospheric emission and use of planetary emission for 
signal calibration. However, because of the poor S/N of the WMAP data, the situation with 
regard to this calibration method is still far from being satisfactory. We also note that an 
attempt to explain the (small) discrepancy between the MITO and OVRO measurements as 
due to confusing CMB anisotropy signal (De Petris et al. 2002) seems less likely in light of 
the WMAP results which imply that such a signal is relatively insignificant. 

In the relativistically accurate calculation of Al, the intensity change in Compton scat- 
tering, its dependence on the gas temperature is no longer linear (Rephaeli 1995b), and can 
be written in the form (Battistelli et al. 2002) 

AT 2k 3 T 3 a T x 4 e x f \kT e , . . v r „, ^ ,i „ 

where T is the CMB temperature, ot is the Thomson cross section, x = hu/kT is the 
nondimensional frequency, n e and T e are the electron density and temperature, v r is the 
line of sight (LOS) component of the cluster (peculiar) velocity in the CMB frame, and 
fi(x) = x(e x + l)/(e x — 1) — 4. The integral is over the LOS through the cluster. Both the 
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thermal and kinematic components of the effect are included in equation (f), separately in 
the first two terms, and jointly in the function R(x,T e ,v r ). Analytic approximations to the 
results of the relativistic calculations have been derived to various degrees of accuracy by 
expansion in powers of kT e /mc 2 and v r /c (e.g., ,, Itoh et al. 1998; Nozawa, Itoh, & Kohyama 
1998; Shimon & Rephaeli 2003; Colafrancesco, Marchegiani, & Palladino 2003). If the gas is 
isothermal and the cluster velocity is sufficiently small so that the velocity-dependent terms 
are negligible (at the relevant frequencies), then AI oc r, and a one parameter fit can be 
performed to the spectral SZ measurements in order to determine the value of r. Since the 
cluster velocity is unknown, the contribution of the kinematic terms can then be treated as 
a source of systematic error. 

The Coma cluster was observed by all X-ray experiments; we use values of the gas 
parameters as deduced by Mohr, Mathiesen, & Evrard (1999) from ROSAT observations. 
The emission-weighted gas temperature is kT e = 8.21 ± 0.16, and the gas core radius, index 
of the /9-model density profile, n e {6) = n e o(l + (6/d c ) 2 )~ 3f3 ^ 2 (where 6 is the angular radial 
variable), and central surface brightness are 9 C = 9'. 97 ± 0'.67, (3 = 0.705 ± 0.046, and 
Sx,o = (4.65±0.14) x 10~ 13 erg s^ 1 cmT 2 arcmin~ 2 , respectively. These values are consistent 
with the results from more recent XMM measurements. 

Best fitting the six data points in Figure 1 by the calculated spectrum (taking kT e = 
8.21 ± 0.16) we obtain the central value of the Thomson optical depth, tq = (5.35 ± 0.67) x 
10~ 3 , where the error reflects a 68% confidence level uncertainty associated with the x 2 mini- 
mization procedure and the uncertainty due to the error in the observed electron temperature 
which contributes ~ 0.10. 



3. The Hubble constant, H 

As is well known, the combination of SZ and X-ray measurements yields a determination 
of the cluster angular diameter distance, d^, from which H can then be deduced by a 
comparison with the theoretical expression for d^- The equation used to determine cIa is 
obtained by combining the formulae for AI and the X-ray surface brightness, Sx, through 
the elimination of the dependence of these quantities on the central electron density. This 
procedure has been discussed in detail in numerous works (e.g., Holzapfel et al. 1997, Hughes 
& Birkinshaw 1998, Furuzawa et al. 1998, Reese et al. 2000, Mason, Myers, & Readhead 
2001) and has by now become standard. Here we employ this method to determine Ho for 
the first time from multispectral measurements of the Coma cluster. 

Assuming an isothermal gas with a /9-model density profile (Cavaliere & Fusco-Femiano 
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1976), the intensity change is AI oc t (1 + (O/Oc) 2 ) 1 ^ 2 ^ 3 ^ 2 ^, and its explicit dependence on 
x and T e is determined by the form of the analytic approximation to the exact relativistic 
calculation. The X-ray surface brightness is 

Sx = A J A +z)4 J ;>> e n,Z|A B ,dC = S Xfl [l + (9/9 c )T^\ (2) 

3 

where z is the cluster redshift, and ( = I /(1a is the cluster path length along the LOS in 
terms of the angular diameter distance. The sum is over the various ionic species (j), with 
densities rij and charges Zj. The thermal bremmstrahlung emissivity in units of erg s _1 cm 3 
is 

A Bj = 1.426 x l(T 27 T e 1/2 f * e~ e G Zj (T e , e)de, (3) 

J € 1 

where e = E/kT e , E is the emitted X-ray photon energy, and G Zj is the Gaunt factor. We 
have used the accurate analytic formula for the Gaunt factor given by Itoh et al. (2000, 
2002). 

A comparison of the deduced value of cIa with the theoretical expression in a cosmological 
model with a cosmological constant then yields the value of H when the values of the 
cosmological density parameters VLm and are specified. The final expression for H is 

H = Anc(l + z) 4 a2T l Xfi6ci { F r n L (z, Q M , n A ), (4) 

where the integrals over the LOS have been performed up to a maximal radius of 10(# c ), by 
defining F n — (1 + ^)"W 2 where £ = C/#c- The function g(z, Q M , f2 A ) is defined as 

f z dz' 

g (z,n M ,n A )= . 3 ==- (5) 

The values Vl M = 0.27 and Vl\ = 0.73 have been adopted in the computation. Doing so we 
obtain 

H = 84 ± 26 km/(sMpc), (6) 

where the error is determined combining the observed SZ (la) uncertainty and the uncer- 
tainties in the X-ray data. These errors contribute comparably to the overall observational 
uncertainty in the deduced value of H . 

To account for the possibility of a nonisothermal profile, we consider a polytropic gas 
model such that the temperature spatial distribution is T e (r) = T e0 [l + (r/r c ) 2 ]~ 3/3 ( 7-1 ^ 2 , 
with the ('polytropic') index 7 considered as a free parameter to be determined (mostly) 
from X-ray measurements. Ideally, when spatially resolved measurements are available of 
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the X-ray spectrum and surface brightness profile, it is possible to determine values of the 
central temperature, central density, core radius, and the indices of the f3 and polytropic 
profiles by direct fits to the data. This complete procedure will be feasible when the full 
XMM data set on Coma becomes available. For now, we use the previously determined 
ROSAT parameters to estimate the implied change in the value of Hq for an assumed range 
of values of 7. For a given value of 7 we determine the central temperature by keeping the 
emissivity-weighted temperature at its observed value (which was deduced by assuming an 
isothermal model). Since the observed surface brightness profile is still fitted by a /3-model, 
but now with the value 0.705 ±0.046 identified as (3 iso (corresponding to the isothermal case, 
7 = 1), it follows that the 7-dependent value of (3 is 

/3( 7 )=4/W(3 + 7)- 

We have considered values of 7 in the range 1 - 5/3. Repeating the above procedure for the 
determination of H for values of 7 in this range, we find that H changes by at most 5% 
with respect to the value computed with 7=1. This variation is well within the observa- 
tional error quoted above, and the overall estimated level of known systematic uncertainties 
stemming from additional simplification in the modeling of the gas, such as spherical sym- 
metry, unclumped gas distribution, negligible impact of CMB anisotropy, and the kinematic 
SZ component, as well as other confusing signals. The additional level of systematic uncer- 
tainties in the value of H (deduced from SZ-X measurements) is estimated to be ~ 30% 
(e.#., Rephaeli 1995a, Holzapfel et al. 1997, Birkinshaw 1999). 

4. Discussion 

Our deduced value for H is about ~ 30% higher than the current mean value from 33 
individual results from the full data set currently available (Carlstrom et al. 2002). Given 
the (relatively) large error, it is fully consistent with the recent value derived by Spergel et 
al. (2003) from the first year WMAP measurements. Clearly, the interest in using the SZ-X 
method to determine H is not diminished by the high quality WMAP result. In addition 
to the need for alternative independent methods to CMB sky maps, measurement of this 
basic parameter at many redshifts and directions on the sky yields important additional 
information on its variation over cosmological time, and test of its predicted isotropy. 

The full potential of the SZ-X method to determine H has not yet been realized. 
As has often been stated by Rephaeli (e.g., Rephaeli 1999), results from this method will 
be optimized when a large sample of nearby clusters will be measured by multi-frequency 
bolometric arrays, and with the availability of high quality X-ray results from the XMM and 
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Chandra satellites. Many ground-based - and stratospheric (such as OLIMPO ; Masi et al. 
2003) - SZ projects will collect sensitive data on a large number of clusters. The overall 
impact of these projects is likely to be larger than a single - even if multiyear - satellite 
project: the main advantages offered by space observations are the very low atmospheric 
background and a more precise calibration based on the modulation of the CMB dipole 
(when compared to calibration by planetary emission). On the other hand, the advantages 
of ground-based observations (from dry sites) are longer observation times, the use of much 
larger telescopes (capable of attaining higher angular resolution), and overall better control 
of systematics that can be obtained with dedicated telescopes. Thus, even though future SZ 
all-sky surveys from space experiments will sample a large quantity of clusters, best results on 
individual clusters are expected from ground-based telescopes. When the sensitivity of space 
observations will reach S/N ~ 100, it will be possible to calibrate ground-based observations 
with the same precision as with space instruments. The uncertainty in determining H will 
then be comparable to that reached by WMAP and optical measurements with the Hubble 
Space Telescope (key project). 

We wish to thank Naoki Itoh for the code of the analytic fitting formula for the Gaunt 
factor. This work has been supported at Rome University by COFIN-MIUR 1998, & 2000, by 
ASI contract BAR, and by a NATO Grant. Work at Tel Aviv University has been supported 
by the Israel Science Foundation. 
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Fig. 1. — SZ spectrum of the Coma cluster. Solid line: best fit spectrum (assuming 
isothermal gas with kT = 8.21 keV) to the combined MITO (diamonds), OVRO (square), 
and WMAP (triangles) measurements, corresponding to r = (5.35 ± 0.67) x 10~ 3 . 
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Table 1. SZE experiment parameters 



Channel 


V 


Au 


f.o.v. (FWHM) 


Sensitivity 


AT 




(GHz) 


(GHz) 


(arcmin) 


(mKs 1 ' 2 ) 


(fiK) 


OVRO 


32.0 


6.5 


7.35 


1.40 


-520±83 


WMAPfy) 


60.8 


13.0 


20 


1.13 


-240±180 


WMAP(W) 


93.5 


19 


13 


1.48 


-340±180 


MITO(l) 


143 


30 


16 


1.21 


-184±39 


MITO(2) 


214 


30 


16 


1.14 


-32±79 


MITO(3) 


272 


32 


16 


0.896 


172±36 



